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Abstract 

Nonlocal coordinate space optical potentials for the scattering of 65 MeV pro- 
tons from nuclei ranging in mass from 6 Li to 238 TJ have been defined by folding 
a complex, medium dependent effective interaction with the density matrix el- 
ements of each target. The effective interaction is based upon solutions of the 
Lippmann-Schwinger and Brueckner-Bethe-Goldstone equations having the 
Paris potential as input. The nuclear structure information required in our 
folding model are the one body density matrix elements for the target and the 
single nucleon bound state wave functions that they weight. For light mass 
nuclei, very large basis shell model calculations have been made to obtain the 
one body density matrix elements. For medium and heavy mass nuclei, a 
very simple shell model prescription has been used. The bound state single 
particle wave functions that complete the nuclear density matrices are either 
Woods-Saxon or harmonic oscillator functions. The former are employed in 
most cases when large basis structure is available. For light nuclei (A < 16) 
Woods-Saxon potential parameters and harmonic oscillator lengths are de- 
termined from fits to electron scattering form factors. For all other nuclei, we 
use harmonic oscillator functions with the oscillator lengths set from an A 1 /® 
mass law. Using this microscopic model, optical potentials result from which 
differential cross sections, analyzing powers and spin rotations are found. In 
general the calculated results compare very well with data when the effective 
interactions are determined from a mapping of nucleon-nucleon g matrices. 
This is not the case when effective interactions built from a mapping of (free) 
t matrices are used. 
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I. INTRODUCTION 



A microscopic, coordinate space model to analyse proton scattering from nuclei has been 
developed |]]|| from that formed earlier by the Hamburg group || . With this model, analyses 
have been made of 200 MeV proton elastic scattering from a large number of nuclei ||, of 
elastic and inelastic p- 12 C scattering over a range of incident energies ||, of charge exchange 
reactions ||, of the structure of neutron rich nuclei such as 6 ' 8 He and 9 ' n Li ||, and, 
very recently, of proton scattering from 3,4 He 0. Those light mass systems usually have 
been considered as few-body problems and the data analysed with few-body methods. We 
note also that microscopic model analyses determined within a momentum space framework 
for elastic proton scattering have been made, with varying degrees of success [|l~0| [3~^1 . The 



energies considered by these microscopic analyses (both in coordinate and momentum space) 
lie in a 'transition' region between low and intermediate energies. For such energies, effects of 
nonlocalities in the effective nucleon-nucleus {NA) interaction must be taken into account. 
Also for these energies, medium dependent effects in the nucleon-nucleon (NN) effective 
interaction upon which that NA interaction is built, are important. When these two facets 
are taken into account in our coordinate space approach, good to excellent agreement with 
data is found to measured elastic and inelastic scattering cross sections, analyzing powers and 
other spin measureables PJTT|. The optical potentials are formed by folding the complex 
effective interaction with nuclear one body density matrix elements (OBDME) and a set 
of single particle (SP) bound state functions. For very light nuclei, modern shell model 
studies |L4| not only specify the OBDME but also which SP wave functions should be used. 
Using this information, calculations of both electron scattering form factors and proton 
elastic scattering observables become predictive; the latter especially when the effective 
interaction and folding procedure are fixed. Otherwise, one can use fits to the longitudinal 
elastic electron scattering form factor to specify the SP wave functions one should use in the 
NA scattering calculations. Those NA calculations then remain predictive as there is no 
adjustable quantity left. For heavier nuclei (A > 16), shell model calculations to date have 
been made only within Ohu space. Therefore we have not used electron scattering form factor 
analyses to select a definitive set of SP wave functions. However, experience suggests || 
that harmonic oscillator (HO) functions with oscillator lengths following a simple mass rule 
(b = A 1 / 6 ) should suffice. 

Whatever the choice of structure input, the folding process leads to NA optical poten- 
tials that are nonlocal because of exchange (antisymmetry) contributions. In our model, 
antisymmetrization of the proton-target system has been taken at the two-body level only, 
i.e. we have used a fully antisymmetrized A-nucleon target wave function and antisym- 
metrized each projectile-target nucleon pair. In the past, the resultant nonlocality of the 
optical potential either was ignored or was approximated by an equivalent local form ||. 
Our calculations of 200 MeV proton-nucleus (pA) scattering gave excellent predictions of 
observables but only when the complete integro-differential forms of Schrodinger equations 
were used M. 

As with our study at 200 MeV, we consider herein only the elastic scattering channel but 
take the spin rotation, R, into account along with the cross sections, da/dQ, and analyzing 
powers, A y . Specifically we have considered 50 targets, namely the Op-shell nuclei 6,7 ' 9,11 Li, 
n B, 12 C, and 16 0; the IsOrf-shell nuclei 20 Ne, 24 Mg, 28 Si, and 32 S; the IsOd proton-0/lp 
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neutron shell nucleus, 40 Ar; the 0/Op-shell 4 o,42,44,48 Ca; 4 6 ,48,50 Ti; 52 Cr> 54,56^ 59 CO; and 
58,60,62,64 Ni . the 2s ld0g neutron shell nuclei 89 Y, 90 Zr, 98 ' 100 Mo, and 118 Sn; the lanthanide 
nuclei, 144 < 152 < 154 Sm, 160 Gd, 164 Dy, 166 > 168 Er and 174 ' 176 Yb; the 0h n/2 proton-0i 13/2 neutron 
shell nuclei 178 > 180 Hf, 182 > 184 W, and 192 0s; the mass 200 nuclei, 208 Pb and 209 Bi; and the 
actinide pair, 232 Th and 238 U. Our predictions of elastic scattering from these diverse mass 
targets are compared with 65 MeV data in all cases with the exception of 9,11 Li; the experi- 
mental data for which come from the elastic scattering of radioactive beams of those lithium 
isotopes from hydrogen. Inverse kinematics makes those experiments equivalent to 60 and 
62 MeV protons scattering from 9 Li and n Li, respectively. 

We compare predictions obtained from the optical potentials defined for each target with 
the proton elastic scattering data that is available, but only to 80° in the centre of mass. In 
general, the cross-section magnitudes by that scattering angle are typically ~ 1 mb/sr. We 
do not expect the approximations needed to make our model practical would be appropriate 
necessarily in making a prediction of the scattering of lesser magnitude. 

The paper is arranged as follows. The procedure for obtaining our microscopic optical 
potentials is outlined in Section [II], as are the amplitudes by which we obtain the proton 
scattering observables. In Section |T| we present and discuss the results for the scattering 
of 65 MeV protons from those nuclei included in the present study. Conclusions are drawn 
in Section [TV]. 

II. THE MICROSCOPIC OPTICAL POTENTIAL 

We develop first the folding procedure by which the nonlocal optical potentials are spec- 
ified. From those, the effective NA interaction is obtained and we define the amplitudes 
that specify the scattering observables. 



A. The folding process 

In a representation with r, r' denoting relative coordinates between a colliding pair of 
particles, the Schrodinger equation describing their scattering by a local Coulomb, Vc(r), 
plus a nonlocal hadronic (optical) potential takes the form 



^V 2 - V c (r) + E 



2/j 



*(r) = J U(r,r')^(r') dr' . 



(1) 



This reduces by using the partial wave expansions, 



im 



(2) 



and 



U(r, r') = £ ^\(a)iX(fir 

Im 



(3) 



to a set of integro-differential equations 



3 



dr 2 



1(1 + 1) 



- V c (r) + E 1 ui(r) = J™ W t (r, r')^(r') dr' . 



(4) 



The Wi(r,r f ) contain both the local interaction and multipoles of the nonlocal interaction. 
Note that for simplicity, we have suppressed all terms due to the intrinsic spin of the system. 
We seek solutions for NA scattering and determine the nonlocal interactions, Una(t, r'), 
at 65 MeV in particular, by folding effective NA interactions with the relevant structure 
information. Thus we obtain the appropriate NA effective interaction for each nucleus in 
our investigation from 6 Li to 238 U. At this particular energy one may anticipate stronger 
influences in analyses of the medium effects defining the effective NN interactions and of the 
nonlocalities in the optical potentials arising from the folding process than in studies of 200 
MeV scattering [[| . With the optical potential in this coordinate space representation, we use 
the program DWBA91 [15| to solve the set of partial wave integro-differential Schrodinger 
equations. That code has the further useful attribute that it can be used to evaluate distorted 
wave born approximated (DWBA) amplitudes for inelastic scattering, given the appropriate 
OBDME and SP wave functions. The microscopic optical potentials are used therein to 
determine the distorted waves and the medium modified effective NN interaction is used 
as the transition operator. Data from inelastic scattering of 200 MeV protons from 6 ' 7 Li 
and 12 C have been well reproduced by using that procedure [p||16! further justifying the 
scattering theory formulated in coordinate space. 

To define the nonlocal interaction for NA scattering in a full folding model, we must 
accommodate antisymmetry between the projectile and every nucleon specified with the 
internal nuclear wave function. We must evaluate multi-particle matrix elements of the 
form 



U pA = (tf 



A 

E Vno 



N=l 



(5) 



with (R| \& (1 . . . A)) being the many-body wave function for the ground state of the target 
and '0' denoting the projectile coordinates. As all nucleons in the target are equivalent, it 
is useful to choose a specific entry ("1") and write 



U pA = A(*(l...A)\V 10 \V(l...A)) 

With the many-body state expanded in cofactors, 

1 



|*(1. ..A)) 



52\(p am (l))a am \V(l...A)) 



where a specifies the set {n, l,j, (}, and ( is the isospin projection, Eq. @ becomes 



pA 



E (* 



(0))} 



(6) 



(7) 



(8) 



ama'm' 



when the required antisymmetry with projectile and struck nucleon is taken into account. 
The nuclear structure information required to evaluate the optical potentials are many-body 
matrix elements of the particle-hole operators. They are defined by 
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mm'MiMf 
Paa'JiJf 



E 



j-m 



(jmj' —m'\I — N) 



i,n y^Jf + 1 
x (J t M t IN\J f M f )S aa/I 

where the OBDME, S aa 'i, are (with d am = (— l) J-m a a _ m ) 



(9) 



Saa'I 



J/ 



(10) 



in the elastic scattering case (from a target with spin J). If J is nonzero, multipoles from 
to I m ax{= 2 J) contribute. Scattering from light odd mass targets show that these must 
be included []16| . As even-even nuclei all have ground state spin-parities (J 7 ") of + , the 
required OBDME simply are the monopoles 



S, 



aa'O 



aL x d c 



(11) 



While the angular momentum selection rules require I = I' and j = j', there is no such 
restriction on the principle quantum number. The cases where n ^ n' signify purely radial 
excitations which can only occur in spaces allowing for cross shell transitions. In those 
instances, the specification of the full density requires the inclusion of those off-diagonal 
elements. The diagonal elements reduce to the shell occupancies as 



a a'm' a ctm 



Qa > mm' O act' > 



where a aa is the fractional shell occupancy as a fully occupied shell has a 
(diagonal) OBDME are given by 



SaaO — Y 2 j + 1 

The optical potential given by this folding process takes the form 



(12) 

1. Thus the 
(13) 



U( ri ,r 2 ;E)= £ 



fl(ri ~r 2 ) / p* a/rn ,(s)U D (R ls ,E)Lp arn (s) ds 
+ cp* alm ,(r 1 )U E %R 12 ,E) lfiam (v 2 ) 



(14) 



where Ru = |ri — ra|, and U and U are appropriate combinations of the multipoles of 
the effective interaction for the direct and exchange contributions to the optical potential 
respectively |3,15|. One then applies the partial wave expansions, Eq. (^). 

A feature in our process of analysis is that all details required to make the calculations 
are preset. In the cases of 6 ' 7 Li |[l(| and of 12 C @, nuclear structure information was taken 
from shell model calculations in complete multi-au; spaces with the SP wave functions set 
by fits to the longitudinal elastic electron scattering form factors. Thus the results for both 
the elastic and inelastic scattering of 200 MeV protons from those nuclei were obtained 
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from single calculations. No adjustments (such as those attributed to core polarizations) 
were needed and most proton excitation data were well described. With the microscopic 
(coordinate space) model for 200 MeV pA scattering established for a set of Op-shell nuclei, 
an extensive analysis for many more nuclei with masses up to 208 Pb was completed In 
almost all cases single calculations led to excellent reproduction of data. Most recently ||, 
wave functions from very large space shell model calculations of 3 He were used in successful 
predictions of the elastic scattering of 200 MeV protons from that few nucleon system. In 
that case the shell model interaction used also fixed the SP wave functions so the analysis 
of the electron scattering form factor was also predictive. 



B. The effective interaction 

The folding procedure to define the optical potential requires the specification of the NN 
effective interaction in ST-channel form and in coordinate space. This effective interaction 
we take as a mixture of central, two-body spin-orbit and tensor force attributes, each having 
a form factor that is a set of four Yukawa functions with complex coefficients, i.e. 

9T ) ^E) = Y.Si ) [p{r\E} e ——. (15) 
k=i ' 

Therein [p(r),E] are complex strengths that vary with projectile energy and nuclear 
density, [i® are the inverse ranges of the Yukawa functions with k the index for those inverse 
ranges. In principle, the number of strengths and inverse ranges (n,) chosen can be as large 
as one likes, though for all operators = 4 seems to be sufficient to reproduce accurately 
the half-off-shell g matrices for laboratory energies between 50 and 400 MeV jl]. Singular 
valued decomposition was used to optimise those inverse ranges and coefficients so that 
double Bessel transforms of the effective interaction map accurately to an appropriate set of 



infinite nuclear matter g matrices |~7|,[18| obtained from solution of the Bethe-Brueckner- 
Goldstone (BBG) equation, 

g { L J L s , T \ P ', P ;k,K,k F ) = vt J J T \p,p') 

+ -E I" V U ST) (P'> ff) [«] 9i J J T \q,P; K K, k F ) q 2 dq, (16) 
n l Jo 

where 

in which Q(q, K, kf) is an angle averaged Pauli operator with an average center of mass 
momentum, K [|17],|18|]. The energies in the propagators of the BBG equations include 
auxiliary potentials, U, and are defined by 



E(q,K,k f ) = [q 2 + K 2 ) + (^) {U (|q + K|) + U (|q-K|)} , (18) 
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wherein U is the first order mass operator term. We take v[jf T ^ (p, p') to be the Paris NN 
interaction |]19||, although there is little difference if one starts with the Bonn NN potential 



0. Details of the techniques involved are given elsewhere fl7jJn 

Given that the NN g matrices are most easily specified in momentum space and the 
effective interaction form is an approximation, it is sensible to seek to analyse NA elastic 
scattering with a momentum space solution of the Schrodinger equation. Such studies have 



been made using credible g matrices as input [0],P!2| and those results reflect reasonable 
agreement with the data. They also confirm the need for inclusion of medium effects for 
low and intermediate energy NA scattering. Of note from the momentum space calcula- 
tions is the observation that off-shell Kowalski-Noyes /-ratios of the t and g matrices vary 
quite similarly |Tj|. The major effect of medium modifications (Pauli blocking and energy 
denominators) are to vary the on-shell values. 

The choice of the Yukawa forms for the effective interaction is extremely advantageous 
when it comes to evaluating the optical potential. Fourier transformation of each of the radial 
components of the effective NN interactions gives a simple multipole form, so that for the 
central terms, the double Bessel transformation that leads to each term in the nonlocal 
interaction can be solved analytically, taking the form 

W?\r u r 2 ) oc (z/i (fe) r>)^ (z^*V<) , (19) 

where r< and r> are the lesser and greater of ri and r 2 respectively. 



C. The observables 



There are diverse observables for the scattering of polarized protons from an unpolarised 
target. While one may define differing sets, we consider that which involves differential 
cross section, da/dQ, analyzing power, A y , and two Wolfenstein spin rotations, A and R. 
These measureables are defined in terms of scattering amplitudes F{9) ||20|| . Since the data 
we investigate were obtained using polarized projectile protons, this amplitude is a 2 x 2 
matrix, 



where 



and 



F{6) 



l 



f- H (0) = fc{e)%_{6) + -L £(2 J + l)e-^ (Sj - Sj) r^\(6). 

2 2 2 2 4:IK ^ V '22 



(20) 



(21) 



(22) 



Here rWi(0) are the rotation matrices, o"j are the Coulomb phase shifts, fc{Q) is the point 

2 ± 2 

Coulomb scattering amplitude, and Sf are the S matrices for J = L ± |. With these, the 
elastic scattering observables are defined as 
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(23) 



A y {B) = P(9) 



2Im< 






2 

+ 


/wo 2 



(24) 



and 



R{9) = KZ{0) = K* K (0) 





2 


/i-j(9) 2 


f H (0) 


2 

+ 





A{0) = K* s (0) = -KW) 



2Re< 


[/ H <Wj-jW} 




2 

+ 





(25) 



(26) 



Commonly Q, a linear combination of the Wolfenstein spin rotations, is measured. It relates 
to the above by 



Q{6) = R{9) sin^ + A(9) cos9. 



(27) 



III. RESULTS 

In the following subsections, we display the results of our calculations of the elastic 
scattering of 65 MeV protons from many target nuclei and place them in comparison with 
experimental data wherever such data exist. The results are subdivided into four sections. 
First, we present the scattering from the light nuclei, 6 Li to 16 0, for which OBDME have been 
obtained mostly from large space shell model calculations and for which SP wave functions 
have been determined from fits to the longitudinal elastic electron scattering form factors. 
Then we present and discuss the results obtained with medium mass nuclei, from 20 Ne to 
64 Ni. For those nuclei the shell model calculations were performed in complete Ohu spaces, 
except for those in the middle of the fp shell where the dimension of the basis becomes 
prohibitively large for matrix diagonalization techniques. In those restriction on the 

number of nucleons in the 0/r orbit was placed on the model. The oscillator length for the 

HO SP wave functions used was set as A 1 / 6 (in units of fm). In the third section, we present 
the results for the scattering from heavy nuclei, and in which the nucleon occupancies are 
those given by a simple packing model. Again we use HO wave functions with an oscillator 
length given by A 1 ^ fm as the SP functions. In the final section, we consider the changing 
structure of the observables across the mass range. 

For each of the cases discussed, we present the results of two calculations. They differ 
only in the specification of the optical potential. The results obtained with the effective 
interaction built from the NN g matrices folded with known density profiles |21| are dis- 
played by the solid curves while those found by using the free NN t matrices are displayed 
by the dashed curves. For simplicity we designate the results as being obtained from cal- 
culations made using ^-folding and i-folding potentials respectively. The data come from 



Refs. [22-33] with specific reference made in the figure captions. 
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A. Light mass nuclei (A < 16) 

The results of our calculations of the elastic scattering of 65 MeV protons from light 
mass nuclei are shown in Figs. [I] and |2|. The differential cross sections, da/dQ, and analyzing 
powers, A y , are displayed in the left and right panels, respectively The target is identified in 
each segment of the diagrams. The ground state wave functions were obtained, for the most 
part, from shell model calculations made within a complete (0 + 2)%uo model space. The 
exceptions were 6 ' 7 Li, for which the wave functions were obtained in a complete (0 + 2 + 4)^c<j 
space [[Tj|. The SP wave functions were assumed to be of Woods-Saxon (WS) form and were 



chosen to reproduce the longitudinal elastic electron scattering form factor data. For 16 0, 
HO functions were used. 

The differential cross sections and analyzing powers for the elastic scattering of protons 
from four lithium isotopes are displayed in Fig. [I]. In the case of 6 Li, our prediction of the 
cross section made with the g-folding potential has the correct shape but it is slightly weaker 



than the data. This is not the case at 200 MeV ||16|| , where the data are well reproduced. 
The slight discrepancy at the lower energy may be a consequence of the deformation of the 
target and the failure of the shell model, even in the (0 + 2 + 4)hu space, to reproduce 



the cluster nature of the Li nucleus Yet, although Li does not have a large central 



density, the changes as a result of the medium modifications contained in the g matrix 
are quite significant. This is observed in comparison to the results obtained using the t— 
folding potential. The measured cross section is not reproduced as well with that result 
as with the ^-folding one, although the discrepancies are not large. However, as noted 
previously f|,|ll|,il], the analyzing powers are very sensitive to the details of the density 
in the folding. We find this to be the case again, and at 65 MeV, with all of the light 
nuclei. Our results of scattering from 7 Li are better, as now our prediction for the cross 
section obtained using the ^-folding potential are only marginally weaker than the data. 
The analyzing power is also well reproduced out to 80°. 

For 9 Li, the medium modification effects on the cross section are quite pronounced. In 
this case, it is essential to use the results from ^-folding. Currently there are no data for 
the spin observables or for form factors from electron scattering for the exotic nuclei so that 
the proton cross section represents the best available test of putative structures. In this 
case, as for n Li, the density profile was obtained directly from the shell model ground state. 
Also, in absence of any electron scattering data, we used WS SP wave functions appropriate 
for electron scattering from 9 Be. With such a specification, the result obtained from the 
^-folding potential is remarkably good. 

The final results shown in Fig. p] are for 11 Li as the target. This nucleus has a halo 
distribution associated with the very loosely bound valence neutron pair. To allow for this 
extension in the density we use WS SP wave functions with the binding energy for the Opi 
and higher shell neutrons set to 500 keV ||. The WS functions assumed for the Osi and 
0»3 orbits were those used for 9 Li, consistent with the model for n Li of a two-neutron halo 

2 

outside a 9 Li core. With that prescription, our predictions adequately reproduce the data, 
although it slightly over-estimates the cross section at the small scattering angles. Even so, 
using such a simple model to describe the halo nature of this nucleus produces results that 
give confidence to use proton scattering to study the matter distributions of other neutron 
and proton rich nuclei. 
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We note in passing that for a target with non-zero spin it is essential to include all 
multipole contributions in the scattering, as is the case with all the lithium isotopes. In 
Eq. fllOD, there are 2 J + 1 possible multipoles in the specification of the ground state. The 
contributions of non-zero multipoles are not large, but they are important. The analyzing 
power is especially sensitive to them, as is most dramatically seen in the results for proton 
scattering from 9 Be P,35|. 

The results for scattering from 11 B, 12 C and 16 are presented in Fig. f^. For 11 B, our 
predictions compare least favourably with the data, although the cross section and analyzing 
power still are well described out to 40°. Above this, the predicted cross section falls at a 
slightly greater rate than that indicated by the data. However, the effects of the medium 
modifications in the effective NN interaction are essential for any adequate description of 
the analyzing power. With the other nuclei, the g-folding model predictions are all very 
good. They are clearly better representations of the data, especially the analyzing powers, 
than are the results found with the t-folding model. Note that the range involves at least 
three orders of magnitude in the cross section. 



B. Medium mass nuclei (16 < A < 64) 

While for light mass nuclei it is now possible to make shell model calculations using 
large and complete spaces, for IsOd and Oflp shell nuclei, such "no-core" determination of 
the shell model wave functions as yet are not feasible. The size of the spaces necessary are 
prohibitively large. For a number of nuclei, however, it is still possible to construct complete 
Ohu) shell models for use in analyses of elastic scattering data. Those models still require 
core polarization corrections and and hence the use of the density matrices in analyses of 
electron scattering data may not be as good a check on the SP wave functions as with the 
Op-shell nuclei. Therefore we have assumed HO functions with A 1//6 fm for all subsequent 
calculations of the proton scattering observables for nuclei with A > 20. 

The results for the elastic scattering of 65 MeV protons from 20 Ne, 24 Mg, 28 Si and 32 S 
are presented in Fig. |3|. A complete Ohu (IsOd) shell model calculation was performed for 
all nuclei. The results of our g-folding model calculations reproduce the cross-section data 
for the scattering from 24 Mg and 28 Si well and at all scattering angles to 80°. For 20 Ne, 
the description is accurate to 60° at least. However for scattering from 32 S, the level of 
agreement with the cross-section data is less than satisfactory. Also in all cases, the level of 
agreement between the results of our calculations using the t-folding model and the data is 
poorer. This poor reproduction we find from our g-folding model to the 32 S cross-section 
data cannot be due simply to the choice of oscillator length. The cross-section data for 32 S 
have magnitudes considerably larger than those of 28 Si (and of 40 Ar that we show next) for 
small scattering angles and vice-versa for larger scattering angles. It is known that these 
IsOd shell nuclei are deformed and that 32 S is different to the others. That is evidenced by 
32 S not having the distinctive splitting of the giant dipole resonances that occur in 24 Mg 
and 28 Si. The different deformation of 32 S might explain the difference we see in the quality 
of reproduction of the scattering data. Certainly when a phenomenological optical model 
analysis was used to fit the same scattering data, the parameters required to fit the 32 S data 



were quite different from those found with data off the neighbouring nuclei [26|. As with 



the results presented for the light nuclei, the differences between the results obtained using 
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the g and t-folding models are far more significant in the analyzing powers. In all cases, 
those results obtained with the g-folding models reproduce the data well. We also note that 
there is a definite trend in the size and shape of the data as one increases the mass of the 
target. The data indicate a sharp rise from at 0° to a maximum near 1 followed by a fall 
and a second peak of similar character. The first peak becomes more forward-peaked as the 
mass of the target increases (60° for 6 Li, 30° for 28 Si). This could be used as an indicator 
for optimum SP wave functions, as small changes in the HO parameter produce shifts in the 
position of this peak in the analyzing power. 

The results for the elastic scattering of 65 MeV protons from 40 Ar and 40 Ca are displayed 
in Fig. The analyzing powers have small differences with the most notable being a small 
shift in the angles of structures. The differential cross sections also are different with the 
position and size of the prominent peak being the primary effect. For 40 Ar that feature occurs 
at 33° and has a value of 60 mb/sr. For 40 Ca, the result is 44 mb/sr at 36°. For both these 
nuclei, a 'packed' model of their structure has been used in our analyses. Specifically 40 Ca 
has been taken as a doubly closed shell nucleus while 40 Ar has been taken as two proton 
holes (in the ob-subshell) with two extra core nucleons (in the fr orbit) on that doubly 
closed shell description. The agreement with data for the scattering from 40 Ar is good to 
60°. The differences between the Ar and Ca results are consistent with observed differences 
in the data sets and reflect the effects of surface contributions revealed by a change of basic 
structure from the closure of a major shell. 

In the case of 40 Ca, our t-folding model results are quite similar to those obtained recently 
in a momentum-space framework [13|, while our ^-folding model results agree quantitatively 



with those found using a g matrix in another momentum space calculation [ll]. These 
differences emphasise that one cannot neglect the importance of the medium in specifying 
the effective interaction, whether it is for scattering at 200 MeV H or at 65 MeV |34| . 

The results of our calculations using the g- and t-folding models are compared with 
data for the case of scattering from the heavier calcium isotopes in Fig. [5[ As with the 
lighter nuclei, we again note that the results found using the optical potentials specified in 
the ^-folding model are in excellent agreement with the data up to 70°. Again those results 
differ markedly from ones obtained using the optical potentials in the t-folding model. Also, 
the trend in the analyzing power noted earlier is observed in the calcium isotopes. However, 
there are more pronounced minima in the cross sections, and the analyzing powers show a 
new minimum developing at small scattering angles. 

The results for the scattering from 46 > 48 > 50 Ti and 52 Cr are displayed in Fig. |6]. The 
agreement between the data for the scattering of 65 MeV protons from the titanium isotopes 
and the results from our g-folding model calculations is quite good to 70°. That level of 
agreement is not observed in the case of scattering from 52 Cr; the data are under-predicted 
above 45°. However, the shape indicated by the data is by far better reproduced by the 
^-folding result than the t-folding one. This under-prediction at large angles is observed 
in the scattering from 54 > 56 Fe and 59 Co, displayed in Fig. [7], and also in the scattering from 
the Ni isotopes, displayed in Fig. |8|. In these complete OHlu shell model calculation 

is not possible with the standard diagonalization techniques; those are only calculable in an 
approximate fashion using Monte Carlo techniques. For the present calculations, the model 
space in the shell model calculation was restricted to close the 0/t orbit. This could cause 
the observed discrepancy with the data. Also, the choice for the oscillator length might not 
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be optimal. One could choose to set b by a best fit to the cross-section and analyzing power 
data, as we did in our 200 MeV analysis [g. Nevertheless, in every case, the results obtained 
in the g-folding model are in much better agreement overall with the data, especially with 
the analyzing powers, with the level of agreement being actually quite good. 

The trends in magnitudes and shape of the data with mass remain and are enhanced with 
this set of nuclei. Notably the first minimum of the cross section becomes more pronounced 
as does the first minimum in the analyzing power. 

C. Heavy nuclei (A > 64) 

For nuclei with mass in excess of 70, shell model calculations, even in restricted spaces, 
generally are not feasible. The increase in the number of valence nucleons, and an increase 
in the number of SP states that may be occupied, cause the dimension of the model space 
to become prohibitively large. Therefore, we have chosen a simple packing model to specify 
the OBDME. In that model, nuclear shells are filled in sequence from the lowest lying shells 
to the Fermi level. As the results of our calculations are not i particularly sensitive to the 
diffuseness in the nuclear density in this mass region, this is not a bad approximation. A 
more important feature is the choice of the oscillator length for the single nucleon bound 
state functions. Again, we have chosen b = A 1 ^ fm for all shells. A more reasonable 
approach would be to assume a different oscillator parameter for the protons and neutrons. 
By that means, the proton and neutron total distributions could be kept similar despite the 
neutron excesses. 

The results for the scattering from 89 Y, 90 Zr and 98 ' 100 Mo are displayed in Fig. [5]. For 
the four cases presented, the cross-section data are well reproduced by our ^-folding model 
predictions to 50°. Thereafter, our results slightly underestimate the data in the region of 
the minimum at ~ 55°. While the t-folding model results give similar shapes for the cross 
sections, the second minimum at 35° is an order of magnitude deeper than that observed 
and also as predicted by the ^-folding model result. The differences between the two models 
are far more striking in the comparisons of the analyzing powers. The ^-folding results are 
in very good agreement with the data to 50°. The results obtained with the t-folding model 
definitely are not. The latter do not reproduce the shape or the magnitude of the data. 
It is interesting to note that the region in which the analyzing power is underestimated by 
the (^-folding results is also that in which the cross section is underestimated. Since the 
analyzing power is scaled by the cross section, an improvement in the level of agreement in 
the cross section in this region may also produce an improvement in the analyzing power. 

Our ^-folding model results are compared with data for the scattering of 65 MeV protons 
from 118 Sn and 144 > 152 > 154 Sm in Fig. pi]. For these cases, the first two minima in the cross- 
section data are very well reproduced by the g-folding model results as is the third minimum 
in the 152 ' 154 Sm data. The same level of agreement is achieved in the analyzing powers 
between the data and results from the ^-folding calculations. The same dramatic difference 
between the g- and t- models is observed in the analyzing power as was the case with the 
results in the mass- 90 region. 

We show comparisons of our model predictions with data for the scattering from nuclei 
ranging between 160 Gd and 180 Hf in Figs. [11] and [12]. The level of agreement in the cross 
section between the data and the ^-folding model results is again very high. The results 
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obtained using the t-folding model have a tendency to underestimate the data, especially 
at the minima above 40°, and do not predict the locations of those minima. Nevertheless 
the shape of the cross-section data generally is reproduced. Those levels of agreement 
are not reflected in the comparisons of analyzing powers. The ^-folding model predictions 
are results that generally reflect the data, although there might be some suggestion from 
the comparisons with the low angle data that larger oscillator lengths are preferable. The 
t-folding model ones reproduce neither the magnitudes nor the shapes of the data. 

The data for the scattering from 182 ' 184 W, and 192 0s are compared with our g— and t- 
folding model predictions in Fig. [T3|. The level of agreement with the cross-section and 
analyzing power data for the scattering from the W isotopes is as observed for the cases 
discussed already. The results for the scattering from 192 Os are a little perplexing. The 
data suggest a somewhat weaker cross section than those for the elastic scattering from 
neighbouring nuclei. Yet the shape and magnitude of the analyzing power is similar. Our 
^-folding results overestimate the cross section by 60 — 70% at low scattering angles. 

We compare the results of our microscopic optical model calculations with the data for 
scattering from nuclei with A > 200 in Fig. [H[ All the data, for both the cross section and 
analyzing power, show similar structure, with which our g-folding model results agree well 
in both shape and magnitude. However, the slight differences in the positions of the minima 
between the cross-section data and our predictions suggest that the choice of oscillator 
wave functions is less than optimal. Yet the comparison with the analyzing power data is 
remarkably good in tracking the shape and positions of the minima. It is clear once more 
that the effects of the nuclear medium in defining the effective NA optical potential are quite 
important. While the t-folding model results track the positions of the maxima and minima 
in the analyzing powers to some extent, they fail to reproduce the observed magnitudes. 
This is also reflected in the cross sections, where the positions of the minima also are not 
reproduced. 



D. Mass dependencies of spin observables 

In Figs. [15] and [16] we display the analyzing powers and spin rotations (R) for the scatter- 
ing of 65 MeV polarized protons from a set of 8 nuclei ranging from 12 C to 208 Pb. The curves 
represent the same model predictions as given in the preceding figures. One can see from 



Fig. [15] that the structure of the measured analyzing power changes in a consistent way as 
the mass of the target increases. Indeed, the depth of the first minimum and the sharp rise 
between this and the next maximum becomes more pronounced with target mass to 118 Sn. 
As the target mass increases, the magnitude of this minimum and of the following maxi- 
mum are greatly suppressed. The absolute value of those minima and maxima approaches 
unity for the heaviest nuclei. As noted in the discussions of individual results, the ^-folding 
calculations predict observed analyzing powers very well and at scattering angles for which 
the cross-section data (usually of magnitude > 1 mb/sr) also are well reproduced. 

There is a limited set of measurements of the spin rotations for the elastic scattering 



of 65 MeV polarized protons from nuclei. These are displayed in Fig. [16]. For the lightest 
four nuclei presented, the comparison between the (/-folding model results and data is quite 
good. Those results from the t-folding model calculations do not match observation as 
adequately. The differences between the g- and t-folding models are far more pronounced 
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with the scattering from the heavier targets. While the ^-folding model results again give 
reasonable agreement with the data, the t-folding model results fail to reproduce both shapes 
and magnitudes. Indeed, that model predicts a maximum at ~ 20° in the spin rotation for 
the scattering from 90 Zr and 118 Sn; the data indicate a minimum at that angle and such is 
predicted by the (/-folding model. While there are no data for the scattering from 152 Sm, 
that difference between the models is also apparent. Also by this mass, the first minimum 
has almost disappeared. In the case of scattering from 208 Pb, the result obtained from 
the ^-folding model reflects the shape of the data, although it overestimates the observed 
magnitude above 30°. A possible improvement to this result is to use a different set of SP 
wave functions, as is discussed below. 



E. Effect of varying the oscillator length 

For the descriptions of the scattering for all nuclei up to and including 20 Ne we have 
used SP wave functions which have been set by fits to electron scattering data. For heavier 
targets, SP wave functions of HO form with b = A 1 / 6 fm have been used. In general, this 
choice has produced very good results in comparison with data, but it is instructive to 
investigate the sensitivity of our calculations to variations in that choice. For this example, 
we compare the results of two (/-folding calculations with scattering data from 58 Ni and 
208 Pb in Fig. [L7[ Therein, the solid line shows the results of our calculations made using the 
standard value. That is 1.97 fm and 2.43 fm for 58 Ni and 208 Pb, respectively. The dashed 
lines display the results obtained when the oscillator length was chosen to give a much better 
fit to the cross-section data. Note that these choices are predicated on the simple structure 
models assumed for both nuclei. Specific shell effects are expected to have some effect on 
the predictions of cross sections. The revised parameters are 1.87 fm and 2.35 fm for 58 Ni 
and 208 Pb, respectively. Those lengths were used for all orbits. The improvement in the 
results with the new SP wave functions is very evident. There is much better agreement 
between the calculation and experiment for both nuclei, although the position of the minima 
in the cross section for the scattering from 58 Ni are greater than observation. Using these 
new sets of wave functions, one observes now much better agreement with the data for both 
spin observables. In the case of 208 Pb, the significant improvement in the spin rotation has 
been achieved with only a 3% decrease in the oscillator length. 



IV. CONCLUSIONS 

Optical potentials for the elastic scattering of 65 MeV protons from nuclei have been 
obtained by folding medium dependent effective NN interactions with a specification of the 
ground state for each nucleus, and also with SP wave functions of either WS or HO form. 
Those optical potentials are complex and nonlocal and the scattering phase shifts and S- 
matrices from which predictions of the measured quantities were obtained, result by solving 
the relevant nonlocal Schrodinger equations. 

We have obtained good to excellent agreement with scattering data from targets ranging 
from 6 Li to 238 U using the optical potentials so defined. The framework by which the 
results were obtained is entirely parameter- free; no adjustment of any part of the input 
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was necessary. Thus the nonrelativistic mean-field theory for pA scattering based upon 
the infinite matter g matrices is reliable for proton energies down to 65 MeV. This gives 
encouragement for these techniques to be used in analyses of data from radioactive beam 
experiments. As the optical potentials are derived directly from the nucleon distributions, 
instead of averaged charge or matter distributions, this would provide detailed information 
on the structures of exotic nuclei, as was the case in the study of 6,8 He and of 9,11 Li ||. 
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FIG. 1. The differential cross sections (left) and analyzing powers (right) from the elastic 
scattering of 65 MeV protons from 6 ' 7 ' 9,11 Li. The data (dots) are compared with the results of our 
microscopic model calculations for the cases when medium effects are included (solid curves) or 
are ignored (dashed curves). Data were taken from Refs. 
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FIG. 2. As for Fig. | but for n B, 12 C and 16 0. Data were taken from Refs. ||-§7|. 
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FIG. 3. As for Fig. but for 20 Ne, 24 Mg, 28 Si and 32 S. Data were taken from Refs. |||27|. 
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FIG. 4. As for Fig. | but for 40 Ar and 40 Ca. Data were taken from Refs. ||-|2|. 
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FIG. 5. As for Fig. |T| but for 42 - 44 ' 48 Ca. Data were taken from Refs. 
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FIG. 6. As for Fig. | but for 4 M8,50 Ti and 52 Cr _ Data were ta k en f rom Re f s . ||J2|. 
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FIG. 7. As for Fig. g] but for 54 > 56 Fe and 59 Co. Data were taken from Ref. [§7|- 
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FIG. 8. As for Fig. £] but for 58,60,62,64 NL Data were taken from Ref 
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FIG. 9. As for Fig. | but for 89 Y, 90 Zr and 98 > 100 Mo. Data were taken from Ref. |7). 
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FIG. 10. As for Fig. @ but for 118 Sn and i44,i52,i54 Sm Data were taken from Refg |||^|||. 
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FIG. 11. As for Fig. Q but for 160 Gd, 164 Dy and !66,l68 Er Data were taken from Refs> ||^|. 
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FIG. 12. As for Fig. | but for 174 > 176 Yb and ^-^Hf. Data were ta k en f rom Re f s . |||| 
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FIG. 13. As for Fig. | but for 182 > 184 W and 192 0s. Data were taken from Refs. |3(|||. 
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FIG. 14. As for Fig. | but for 208 Pb, 209 Bi, 232 Th and 238 U. Data were taken from Refs. 
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FIG. 15. The analyzing power, Ay, from the elastic scattering of 65 MeV protons from 12 C, 
16 0, 40 Ca, 58 Ni, 90 Zr, 118 Sn, 154 Sm and 208 Pb. The data (dots) are compared with the results of 
our microscopic model calculations for the cases when medium effects are included (solid curves) 
or are ignored (dashed curves). Data were taken from Ref. 
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FIG. 17. The differential cross sections (top), analyzing powers (middle) and spin rotation 
(bottom) from the elastic scattering of 65 MeV protons from 58 Ni and 208 Pb. The data (dots) are 
compared with the results of our microscopic model calculations when differing oscillator lengths 
for the bound state wave functions are used. The solid lines represent the b = A 1 ^ choice displayed 
previously (1.97 fm for 58Ni, 2.43 fm for 208 Pb), while the dashed curves represent fitted values 
for b (1.87 fm for 58Ni, 2.35 fm for 208 Pb). Data were taken from Refs. (2|j27]. 
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